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ABSTRACT 

We have studied the abundance of oxygen in the IGM by analyzing VI, C IV, Si IV, and H I 
pixel optical depths derived from a set of high-quality VLT and Keck spectra of 17 QSOs at 2.1 < 
z < 3.6. Comparing ratios tovi/'?~Civ('?civ) to those in realistic, synthetic spectra drawn from a 
hydrodynamical simulation and comparing to existing constraints on [Si/C] places strong constraints 
on the ultraviolet background (UVB) model using weak priors on allowed values of [Si/O]: for example, 
a quasar-only background yields [Si/O] ps 1.4, highly inconsistent with the [Si/O] ps expected from 
nucleosynthetic yields and with observations of metal-poor stars. Assuming a fiducial quasar -(-galaxy 
UVB consistent with these constraints yields a primary result that [O/C] = 0.66 ± 0.06 ± 0.2; 
this result is sensitive to gas with overdensity 5 > 2. Consistent results are obtained by similarly 
comparing tovi/thi(thi) and Tovi/''~Siiv('?Siiv) to simulation values, and also by directly ionization- 
correcting tovt/thi a s function of thi into [O/H] as a function of density. Subdividing the sample 
reveals no evidence for evolution, but low- and high-rui samples are inconsistent, suggesting either 
density-dependence of [O/C] or - more likely - prevalence of collisionally- ionized gas at high density. 

Subject headings: cosmology: miscellaneous — galaxies: formation — intergalactic medium — quasars: 
absorption lines 



1. INTRODUCTION 

The enrichment of the intergalactic medium (IGM) 
with heavy elements has, over the past decade, become a 
key tool in understanding star and galaxy formation by 
providing a fossil record of metal formation and galactic 
feedback. 

Absorption line spectroscopy has revealed, among 
other findings, that the low-density (6 = pj (p) < 10) 
intergalactic medium (IGM), as probed by the Lya for- 
est and through C III, C IV, Si III, Si IV, O V, O VI, and 
other transitions, is at least partly enriched at all red- 
shifts and densities probed. In particular, recent studies 
indicate that: 

• When smoothed over large (~ 10 — 10 2 kpc) scales, 
the abundance of carbon decreases as gas overden- 
sity S does and has a scatter of ~ 1 dex at fixed 
density. There is carbon in at least some gas at all 
densities down to at least the mean cosmic density, 
with the median carbon metallicity obeying [C/H] 
ps -3.5 + 0.65(log(5 - 0.5) at z ps 3 (Schaye et al. 
2003, hereafter Paper II). 

• On smaller (< 1 kpc) scales the distribution of 

1 Based on public data obtained from the ESP archive of obser- 
vations from the UVES spectrograph at the VLT, Paranal, Chile 
and on data obtained at the W. M. Keck Observatory, which is 
operated as a scientific partnership among the California Institute 
of Technology, the University of California, and the National Aero- 
nautics and Space Administration. The W. M. Keck Observatory 
was made possible by the generous financial support of the W. M. 
Keck Foundation. 

2 Department of Physics, University of California at Santa Cruz, 
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metals is less well known, but observations sug- 
gest that the metals may be concentrated in small, 
high- metallicity patches (Schaye et al. 2007). 

• There is no evidence for metallicity evolution from 
redshift z ps 4 to z ps 2 (Paper II) and metals exist 
at some level at z ps 5 — 6 (Songaila 2001; Pettini 
et al. 2003; Ryan- Weber et al. 2006; Simcoe 2006). 

In connection with this observed widespread distribu- 
tion of metals, a general picture has emerged that galac- 
tic winds - driven largely from young and/or star burst 
galaxies - have enriched the IGM. The same feedback 
may account for the dearth of low-luminosity galaxies 
relative to the halo mass function (e.g., Kauffmann et 
al. 1993; Somerville & Primack 1999; Springel & Hern- 
quist 2003), and also for the mass-metallicity relation- 
ship of galaxies (e.g., Tremonti et al. 2004; Erb et al. 
2006). However, a detailed understanding of the vari- 
ous feedback processes is lacking and there are still open 
questions and controversies concerning the time and rela- 
tive importance of the various enrichment processes, and 
concerning the implications for galaxy formation. 

Both theoretical modeling and observations of inter- 
galactic (IG) enrichment are now advancing to the point 
where comparison between the two can provide crucial 
insight into these issues, but this comparison is not with- 
out problems. Two key difficulties concern the ionization 
correction required to convert observed ionic abundances 
into elemental abundances. First, while the oft-studied 
ions C IV and Si IV are observationally convenient, they 
are poor probes of hot (> 10 5 K) gas, because the ion 
factions C IV/C and Si IV/Si both fall dramatically with 
temperature. Thus, the hot remnants of fast outflows 
might be largely invisible in these ions. Second, the dom- 
inant uncertainty in both the absolute and relative abun- 
dance inferences stems from uncertainty in the spectral 
shape of the ultraviolet ionizing background radiation 
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(UVB). 

Analysis of oxygen, as probed by O VI, has the poten- 
tial to shed light on both problems: this ionization state 
becomes prevalent in some of the very phases in which 
C IV and Si IV become rare, and its abundance de- 
pends on the UVB shape differently than those of other 
ions, helping break the degeneracy between abundances 
and UVB shape. The challenge posed by O VI is that 
at z > 2 it is strongly contaminated by both Lya and 
Ly/3 lines, making its identification and quantification 
difficult. Previous studies of high- z oxygen enrichment 
using line fitting (Carswell et al. 2002; Bergeron et al. 
2002; Simcoe et al. 2004) or pixel statistics (Schaye et 
al. 2000a; Telfer et al. 2002; Pieri & Haehnelt 2004) have 
reliably detected oxygen in the IGM, and quantified its 
abundance in relatively dense gas, but have not assessed 
the oxygen abundance with a very large data sample, at 
very low-densities, or in a unified treatment with other 
available ions. 

Here we extend to O VI our application of the "pixel 
optical depth" technique (e.g., Aguirre, Schaye & The- 
uns, hereafter Paper I) to a large set of high qual- 
ity VLT/UVES and Keck/HIRES spectra. The re- 
sults, when combined with previous studies of C IV and 
C III (Paper II) and of Si IV and Si III (Aguirre et. al. 
2004; hereafter Paper III), give a comprehensive obser- 
vational assessment of IG enrichment by carbon, silicon 
and oxygen, with significantly reduced uncertainties due 
to the UVB shape, as well as new data on the importance 
of hot, collisionally ionized gas. 

We have organized this paper as follows. In §§2 and 3.1 
we briefly describe our sample of QSO spectra. The 
analysis method is described briefly in §3.1 and then in 
greater depth in the remainder of §3, with heavy refer- 
ence to Papers I, II and III. The basic results are given 
in §4 and discussed in §5. Finally we conclude in §6. 

All abundances are given by number relative to hy- 
drogen, and solar abundance are taken to be (O/H) = 
-3.13, (C/H) = -3.45, and (Si/H) Q = -4.45 (Anders 
& Grevesse 1989). 

2. OBSERVATIONS 

We analyze 17 of the 19 high-quality (6.6 kms" 1 veloc- 
ity resolution, S/N > 40) absorption spectra of quasars 
used in Papers II and III. The two highest-redshift spec- 
tra used in those previous studies were excluded here be- 
cause the severe contamination of the O VI region by H I 
lines makes detection of O VI nearly impossible and also 
introduces very large continuum fitting errors in the O VI 
region. Fourteen spectra were taken with the UV- Visual 
Echelle Spectrograph (UVES, D'Odorico et al. 2000) on 
the Very Large Telescope and three were taken with the 
High Resolution Echelle Spectrograph (HIRES, Vogt et 
al. 1994) on the Keck telescope. For convenience, the 
observed QSOs are listed in Table 1. 

Regions within Av = 

max(4000, 8 Mpc H(z)/h) kms -1 from the quasars, 
where H(z) is the Hubble parameter at redshift z extrap- 
olated from its present value (Hq = lOOh kms -1 Mpc -1 ) 
assuming (fi m ,fiA) = (0.3,0.7), were excluded to avoid 
proximity effects. Regions thought to be contaminated 
by absorption features that are not present in our 
simulated spectra (e.g., damped Lya systems) were also 
excluded from the analysis. 
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References. — (1) Kim et al. 2002; (2) Kim et al. 2004; (3) 
Kim, Cristiani, & D'Odorico 2001; (4) Rauch et al. 1997. 



Note. — Columns 1 and 2 contain the quasar name and cor- 
responding Lya emission redshift. Columns 3 and 4 contain the 
minimum and maximum absorption rcdshifts considered and Col- 
umn 5 contains the corresponding minimum observed wavelength. 
The last column contains an estimate of the percentage uncertainty 
in the continuum fitting in the O VI region. 



Lyman continuum contamination increases signifi- 
cantly towards lower wavelengths, whereas (as described 
below) our correction for this contamination assumes 
that it is non-evolving. To mitigate this effect, only the 
red portion [> med(z)] of the QSO spectra used in Pa- 
pers II and III is analyzed in this Paper. As in Schaye 
et al. (2000a), this was found to result in smaller errors 
than using the full region. 

Further details concerning the sample and data reduc- 
tion are given in Paper II (§2). 

3. METHOD 

The pixel optical depth method we use for measuring 
O VI is similar to that described in Papers I, II, and III. 
Section 3.1 contains a brief outline of the method; §3.2 
and §3.3 describe continuum fitting and contamination 
corrections, which have been changed slightly from the 
methods described in Papers II and III; §3.4 describes 
tests of the recovery, and §3.5 discusses the ionization 
balance of the relevant species, and describes how ion- 
ization corrections are performed. 

3.1. Overview 

The basic method for analysis of each QSO spectrum 
is as follows: 

1. Optical depths due to H I Lya (A1216 A) absorp- 
tion are recovered for all pixels in the Lya forest 
region, using higher-order Lyman lines to estimate 
optical depths for saturated pixels. 

2. The pixel optical depth at the corresponding wave- 
lengths of the metal lines O VI (AA1032, 1038), 
C IV (AA1548, 1551), and Si IV (AA1394, 1403) 
are recovered, making several corrections to reduce 
contamination and noise. 
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3. The recovered optical depth in one transition is 
compared with that of another, by binning the pix- 
els in terms of the optical depth of H I, C IV, or 
Si IV, and plotting the median (or some other per- 
centile of) optical depth of O VI. A correlation then 
indicates a detection of O VI absorption. An ex- 
ample is shown in Fig. 1. 

As was done in Papers II and III, an identical analy- 
sis is applied to synthetic spectra generated using a cos- 
mological, hydrodynamical simulation, kindly provided 
by Tom Thcuns. For each observed quasar we generate 
50 corresponding simulated spectra with the same noise 
properties, wavelength coverage, instrumental broaden- 
ing, and pixel size as the observed spectra. For each UVB 
model (of which several are used; see below) the carbon 
distribution as measured in Paper II, and the value of 
[Si/C] from Paper III, are imposed on the fiducial spec- 
tra. An oxygen abundance is assigned by assuming a 
constant, uniform value of [O/C]. Ionization balances are 
calculated using CLOUDY° (version 94; see Ferland et al. 
1998 and Ferland 2000 for details). A direct comparison 
of the results from these simulated and observed spectra 
allows for inferences about the distribution of oxygen, 
carbon, and silicon. The same simulation was used in 
Papers I (§3), II, and III, to which the reader is referred 
for details. 

This study employs the identical UVB models used 
in Papers II (§4.2) and III, excluding model "QGS3.2" . 
All models are from Haardt & Madau (2001, hereafter 
HM01) 6 . These have been renormalized (by a redshift- 
dependent factor) such that the simulated spectra match 
the observed evolution of the mean H I Lya absorption 
(Paper II). The fiducial model, "QG", includes contri- 
butions from both galaxies (with a 10% escape fraction 
for ionizing photons) and quasars; "Q" includes only 
quasars; "QGS" is an artificially softened version of QG: 
its flux has been reduced by a factor of ten above 4 Ryd. 
The UVB used in the simulation only affects the IGM 
temperature, and was chosen to match the measurements 
by Schaye et al. (2000b). 

3.2. Continuum fitting 

A major source of error in O VI optical depths is con- 
tinuum fitting in the O VI absorption region, where con- 
tamination by Lya and Ly/3 lines is heavy. To make this 
fitting as accurate as possible, and to furnish an esti- 
mate of the continuum fitting error, we have applied the 
following procedure to the region analyzed for O VI ab- 
sorption (in the case of observed spectra, this was done 
after the spectra had been continuum fitted by eye as 
described in Paper II §2): 

1. The spectral region is divided into 20 A (rest- 
frame) segments. 

2. In segments with large unabsorbed regions, an au- 
tomatic continuum fitting algorithm is applied (see 
§5.1 of Paper II), in which pixels > la below the 
continuum are iteratively removed. 

5 See http://www.pa.uky.edu/~gary/cloudy. 

6 The data and a description of the input parameters can be 
found at http://pitto.mib.infn.it/~haardt/refmodel.html. 



3. In segments without large unabsorbed regions, we 
identify small unabsorbed regions or regions ab- 
sorbed only in Ly/3; the latter are identified by 
superimposing the region of the spectrum corre- 
sponding to Lya absorption. The continuum level 
of the segment is fit by minimizing the deviation of 
identified unabsorbed regions from unit flux, and 
deviation of the Ly/3 regions from the correspond- 
ing scaled Lya features. 

4. A spline is interpolated between the fits to all seg- 
ments and the spectrum is rescaled by this spline. 

This procedure was applied to all observed spectra as well 
as to to one simulated spectrum per observed spectrum, 
where a 10-20% error in the continuum was introduced 
on scales of 1, 4 and 16 segments. The median absolute 
errors remaining after blindly fitting the continua of the 
simulated spectra are given in Table 1 as an estimate of 
continuum fitting errors in the corresponding observed 
spectra. Because the procedure is not fully automatic, 
we were unable to apply it to all of the simulated spectra. 

The region redwards of Lya was fit in both simulated 
and observed spectra using the procedure described in 
§5.1 step I of Paper II. Continuum fitting errors are much 
smaller for this region (~0.01%). 

3.3. Correcting for contamination 

After continuum fitting the spectra, H I (Lya) optical 
depths thi are derived for each pixel between the quasar's 
Lya and Ly/3 emission wavelengths, save for regions close 
to the quasar to avoid proximity effects (see §2). If Lya 
is saturated (i.e., F(X) < 3cr(A), where F and a are the 
flux and noise arrays, see Paper I, §4.1; Paper II, §5.1, 
step 2), higher-order Lyman lines are used to estimate 
Thi- 

Corresponding O VI, Si IV, and C IV optical depths 
(''ovi, T Siiv, i~civ) are subsequently derived for each H I 
pixel. We exclude regions of the quasar spectrum that 
are contaminated by absorption features that are not in- 
cluded in our simulated spectra, such as Lya lines with 
damping wings. For Ts;rv and tqiv, corrections are made 
for self-contamination and contamination by other metal 
lines, as described in Paper I, §4.2. 

As shown in Fig. 1, when plotting each percentile in 
tovi absorption against absorption in some other ion, 
the correlation disappears below some O VI optical depth 
Tmin (corresponding to a value r c in the other ion) that is 
determined by noise, continuum fitting errors, and con- 
tamination by other lines. These effects may then be 
corrected for by subtracting T m ; n from the binned optical 
depths, thus converting most points below t c into upper 
limits. For each realization and for each percentile, we 
compute T m j n as the given percentile of optical depth for 
the set of pixels with optical depth < r c . We use values 
r c = 0.01 when binning in C IV or Si IV, and of t c = 0.1 
when binning in H I, as we never see a correlation ex- 
tending below these values. 7 

The error on T m ; n for an individual realization is com- 
puted by dividing the spectrum into 5 A segments, then 
bootstrap-resampling the spectrum by choosing these 

7 In Paper II wc used functional fits to determine r c . For OVI 
the correlations are generally less strong than for CIV and wc fix 
t c "by hand." 
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Fig. 1. — O VI optical depths as a function of tqiv (left), thi (middle), and rsiiv (right)- From bottom to top, the points represent the 
31st, median, 69th, and 84th percentiles. Data points are plotted with 1 a error bars. The red lines are predictions from simulations with 
<[C/H]) = -3.8 + 0.655, a = 0.70, [Si/C]= 0.77, [0/C]= 0.64. The observed and simulated 31st, 69th and 84th percentiles have been been 
offset vertically by -0.5, 0.5 and 1.0 dex, respectively. 



chunks with replacement, and finally computing the vari- 
ance of T m in as computed from 100 such resampled spec- 
tra. When the realizations are combined, r m i n is instead 
computed as the median among the realizations, and the 
error on this value is computed by bootstrap-resampling 
the realizations. For further details see Paper II, §5.1, 
step 4 and Paper III, §3.4. 

As noted above, Lya and higher Lyman transitions 
heavily contaminate the O VI absorption regime. This 
can add substantial error in the recovered rovi- Two cor- 
rections are made to minimize this contamination. First, 
after recovering thi and rovi i an initial correction is 
made for contamination by higher-order H I lines by sub- 
tracting 

(A) (A / )(/Ly l A L yi//Ly Q ALya), (1) 

where / is the transition oscillator strength, A' = 
(ALya/AL y i)A is the redshifted Lyman a wavelength cor- 
responding to hyi absorption observed at wavelength A 
and i corresponds to first five higher order Lyman lines 
(i.e., Ly/3 A1025 through LyC A930). 

A second correction is made by taking the minimum 
of the O VI doublet, 



tqvi = mm tqvi 



ZoVIaAoVIaT OVIb 
ZoVIbAoVIb 



(2) 



where 'a' and 'b' denote the stronger and weaker dou- 
blet components, respectively. These corrections are de- 
scribed in detail in Paper I, §4.2. 

Another potential contamination issue is that due to 
strong Lya or Ly/3 absorption, some higher percentiles in 
O VI absorption can become dominated by saturated pix- 
els, so that the particular value of the percentile is deter- 
mined by the contaminating lines rather than by the O VI 
distribution. To remove these unreliable percentiles from 
consideration, the average noise a sa t is calculated for sat- 
urated O VI pixels and is converted into a maximum op- 
tical depth (r sa t = — In 3a sat ). After O VI optical depths 



are binned, those percentile bins with rovi > r sat are 
excluded from the analysis. 

3.4. Testing O VI Recovery 

Because there is substantial processing of the recov- 
ered O VI optical depths, it is important to test how 
efficiently the true O VI optical depths are recovered by 
our procedures. To do so, simulations were produced 
just as described in § 3.1, but with (effectively) perfect 
resolution, no noise, and only O VI A1032 absorption 
(and hence no contamination or self-contamination). In 
Fig. 2 the recovered O VI pixel optical depths are plotted 
against these "true" optical depths for a set of 60 simu- 
lated spectra, for two representative QSOs (see Paper I 
for more such tests). Of particular note is the efficacy of 
subtracting the "flat level" T m ; n as described in the pre- 
ceding section. Ideally, T min would be determined using 
pixels with negligible O VI absorption; this is possible 
in the present case (as the true O VI optical depths are 
known) and corresponding results are shown in the left 
panels. In a realistic case, a proxy for O VI must be used; 
in the right panels of Fig. 2, H I is employed and r m ; n 
is computed using all pixels with logrni < — 1. In both 
cases, the subtracted r m ; n values are shown as horizontal 
dashes on the right axes of Fig. 2. 

Overall, we find that using thi to calculate T m ; n is effec- 
tive at recovering rovi > 10 -2 ' 5 for the 31st and median 
percentiles; for higher percentiles the recovery is accu- 
rate only at higher rovi: but the large 'scatter' indicates 
that this is random, rather than systematic error. 

3.5. Ionization corrections 

In Papers I and II it was shown from simulations that 
there exists a tight correlation between thi and the ab- 
sorbing gas density and temperature, which could be 
used to predict an ionization correction (i.e., the ratios of 
O VI/O and H I/H) as a function of density; for details 
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Fig. 2. — The accuracy of recovering the "true" OVI pixel opti- 
cal depth from simulated spectra of Q0109-3518 and Q1422+230. 
For every panel, the recovered OVI optical depth is plotted against 
the "true" optical depth, rtruc- Each panel shows the binned per- 
centiles after T m i n has been subtracted from the recovered rovi 
pertaining to that percentile. The left and right panels calculate 
T min using the log rtruc < — 3 and logrjji < — 1 respectively. These 
are shown as horizontal dashes on the right axis of each plot. For 
log-rovi > — 2 the median OVI optical depths are effectively re- 
covered. 
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see § 6 in Paper I and § 5.1 in Paper II. As noted in Pa- 
pers II and III, this works well for C IV and less well for 
Si IV, due to their mild and strong ionization correction 
dependence on thi, respectively. 

In the upper panel of Figure 3 we provide a con- 
tour plot of the logarithm of the predicted fraction of 
O VI ions versus temperature and density. The middle 
and lower panels show logTovi/ T civ f° r [O/C]=0, and 
logTovi/ r Hi for [O/H]=0, respectively. For photoion- 
ized gas (T < 10 5 K) the O VI fraction is highest for 
— 5 < logrifj < —4, and only weakly dependent on the 
temperature. 

However, for higher densities (and T < 10 5 K), the 
O VI fraction falls quickly, resulting in a very large ion- 
ization correction for H I saturated pixels. At the same 
time, the O VI fraction at high density increases with T 
for T > 10 5 K, so that collisionally ionized gas might be 
detected relatively easily. Therefore at high thi, colli- 
sionally ionized O VI gas can easily swamp photoionized 
O VI. This can be seen most clearly in the bottom panel, 
which shows that for fixed [O/H] and density, tovt/tji 
increases quickly at T > 10 5 K, particularly if the den- 
sity is high. Because our thi — P relation is dominated 
by photoionized gas (with T < 10 5 ), the effects of col- 
lisionally ionized gas are potentially important and are 
discussed at length in § 5.4 below. 

The strategy employed here is to employ our fiducial 
ionization corrections as in Papers II and III, but to rec- 
ognize that at high-density, the results may be signifi- 
cantly affected by collisionally ionized gas. It is impor- 
tant to note that the importance of collisionally ionized 
gas may be underestimated by our simulation, because 
we did not include a mechanism for generating galac- 
tic winds, which could shock-heat the gas surrounding 
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Fig. 3. — OVI ion fraction Nqyi/Nq (top), log tovi/ t civ f° r 
[O/C]=0 (middle), and log rovi /tii for [O/H]=0 (bottom) as func- 
tions of temperature and the hydrogen number density. Solid 
(dashed) contours are for the z = 3 UV background model QG 
(Q). A measurable quantity of OVI could be collisionally ionized 
if T > 10 5 K, which would be misinterpreted as photoionized gas 
during the ionization correction process. This gas would poten- 
tially contain very little associated CIV and HI. 



galaxies. Our simulation does, however, include heating 
by gravitational accretion shocks. 

Once the ionization correction has been determined, 
and corrected O VI optical depths and recovered H I 
optical depths thi have been obtained, the oxygen abun- 
dance can be calculated: 

[O/H] = log f ™ _22_5E^ _ (0/H) (3) 

V t hi (/ Ajovi n ovi n H J 

where /, and Xi are the oscillator strength and rest 
wavelength of transition i, respectively (/ovi = 0.1329, 
/hi = 0.4164, Aovi = 1031.9270 A, A H i = 1215.6701 A), 
and we use the solar abundance (0/H) Q = —3.13 (num- 
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Fig. 4. — The median oxygen abundance as a function of the 
overdensity (bottom axis) or thi (top axis), by applying the ion- 
ization correction of §3.5 and Paper II. For reference, the dashed 
line is [O/H] vs. density, assuming the carbon distribution fit for 
Q1422+230 from Paper II and constant [O/C]=0.64. 



ber density relative to hydrogen; Anders & Grevesse 
1989). An example of the results from this analysis ap- 
plied to the observed spectrum of Q1422+230 is shown 
in Figure 4. 

In Figure 5 we show a test in which we have generated 
simulated spectra using the "QG" ionizing background, 
recovered optical depths, and applied the just-described 
ionization correction to recover the oxygen abundance. 
The true metallicity is given by the carbon distribu- 
tion of Paper II (for the QG background), with a fixed 
[0/C]=log4.5 w 0.65 (i.e., [O/H]=-3.15+0.65£), and is 
shown on the plot as a dashed line. For log <5 < 1.5, 
the O VI abundance recovery is promising: it overesti- 
mates by less than 0.3 dex, and the dependence of 6 is 
reproduced. However, the overestimation appears to in- 
crease for log S > 1.5, reaching approximately 1 dex for 
the highest overdensity bin. The difference in the high- 
versus \ow-S gas can probably be attributed to the colli- 
sionally ionized gas residing in and around dense regions, 
due to gravitational accretion shocks. We should thus 
keep in mind that we expect to overestimate the oxygen 
abundance associated with strong H I absorbers. 

4. RESULTS 

Our basic result will be an estimate of [O/C] for the 
low-density IGM as a whole, computed using four differ- 
ent but consistent methods, which we describe in turn. 

4.1. ro vi versus tciv f or the full sample 

To extract as much information as possible from our 
data we have, as in papers II and III, combined the data 
points obtained from our entire sample. Figure 6 shows 
log tovi / T civ versus logrciv, in bins of z. To generate 
these points, we begin with tovi values binned in tciv 
for each QSO, as in Fig. 1 for Q1422+230. We then 
subtract from each the "flat level" r m i n for that QSO 
to adjust for noise, contamination, etc. (see § 3.1), then 
divide by the central value of the tciv bin. These points, 
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Fig. 5. — Test of the OVI ionization correction from the syn- 
thetic spectra using the full QSO sample (using 50 realizations of 
each synthetic spectrum so as to focus on systematic effects). The 
input line represents the OVI distribution imposed on the simula- 
tions, and the output line is that recovered from the analysis and 
ionization correction, fit to the data at log<5 <1.5; in this density 
range, the recovery is accurate to within 0.3 dex. 



gathered from all QSOs, are rebinned by determining, for 
each t civ bin in Fig. 6, the best constant-level x 2 ~nt to 
all of the points in the specified rcdshift bin. The errors 
represent 1- and 2-er confidence intervals (A^ 2 = 2 and 
Ax 2 = 4) on this fit. 

The plotted lines indicate corresponding optical depths 
from synthetic spectra drawn from the simulation, using 
several UVB models, the corresponding [C/H] distribu- 
tions as determined in Paper II, and a constant [O/C] 
value determined as follows. For each background, we 
generate simulated tovi/ t civ points in the same way as 
we did for the observations, but averaging over 50 simu- 
lated realizations as described in § 3.1. We then calculate 
a x 2 between all valid observed original (not rebinned) 
points and the corresponding simulated points. 8 Because 
we use 50 simulated realizations, the simulation errors 
are almost always negligible compared to the observed 
errors, but they are still taken into account by calculat- 
ing the total x 2 using the formula: 



E 



X, 



obs 



■^sim 



O'obs 



X, 



obs 



where X = tovi/tciv and a is the error in this quan- 
tity. We then add a constant offset to the simulated 
points (which corresponds to scaling [O/C]) such that 
X 2 is minimized. In each panel the lines connect the 
scaled, rebinned simulation points. 

The first evident result is that log rovi /tciv ~ and 
appears to be at most weakly dependent on log tciv i 
from logrciv — 1.5 to 0. This is unlike rsnv, which 
increases by «2 dex in this tciv regime (Paper II). At 
the lowest densities the data exhibits a decline. Com- 
paring these panels suggest there is little dependence on 
redshift in this interval. This can be seen more clearly 
in Figure 7, which show log rovi / Tciv versus z in bins of 
rciv : there is no evidence, in either the simulated or ob- 

8 Even using 50 realizations, it may occasionally happen that 
a simulated bin fails to have enough pixels for at least five real- 
izations, and so is undefined; in this case the observed point is 
discarded as well. 
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Fig. 6. — Rebinned median log tovi / T CIV vs. logrciv in bins 
of z for the combined QSO sample. The first two panels show 
bins centered at z = 2 and 3 with width Az = 1; the bottom panel 
shows combined data for all redshifts. Data points are plotted with 
1 and 2<r error bars. The lines represent corresponding simulation 
points (with errors suppressed and with [O/C] chosen to minimize 
the x 2 ) using different UVB models, as indicated in the legend of 
the top panel. 



served points, for evolution in tovi/tciv for l-5< z <3.5. 

The observed trends in log(Tovi/ r civ) are reproduced 
well by the simulations. Therefore, because tovi/tciv 
scales with [O/C], the offset in tovi/tciv obtained by 
minimizing the % (Eq. 4) against the observations can 
be used to reliably compute the best fit [O/C]. As an 
example, for our fiducial UVB model QG, the simulated 
spectra were generated with [O/C] =0.65 and best fit by 
an offset of +0.01 dex (implying a best-fit [O/C]=0.66), 
with x 2 /d.o.f.=42.7/82. As we found in Paper II and for 
Q 1422+230 above, the reduced % 2 is somewhat low. This 
is due to a slight overestimate of the errors at low-rciv 
(Paper II) and to the fact that the data points are not 
completely independent (because single absorbers con- 
tribute to multiple data points). 

The fitted [O/C] values and corresponding x 2 /d.o.f., 
are listed in Table 2, with errors computed by bootstrap- 
resampling the quasars used in the x 2 minimization. For 
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Fig. 7. — Rebinned median log(rovi/'?"Civ) vs - z m cuts of tciv 
for the combined QSO sample. Data points are plotted with 1 and 
2<x error bars, where green lines denote lower error of -co. The first 
three panels show bins centered at log tciv = —1.25,-0.75 and 
-0.25 with width 0.5, dex; the bottom-right panel shows data for 
all tciv combined. The lines are the same as in Fig. 6. 



our fiducial model, QG, the best fit [0/C]= 0-66±g;g|. 
The quasar-only background Q (which is probably too 
hard; see Paper II) gives a much lower value of [0/C] = 
0.06tQ Qg. The softer QGS backgrounds gives implau- 



sibly high values of [0/C]= 2.231^. Coupling this 
with results from Paper II suggesting that the QGS back- 
ground is unrealistically soft, strongly disfavors this UVB 
model. 

We may also subdivide our sample by redshift to 
test the dependence of [O/C] on these. First, com- 
puting [O/C] using only spectra that have a median 
absorption redshift med(z) > 2.5 (see Table 1) yields 
[O/C]=0.71 ± 0.07, versus [O/C]=0.58 ± 0.10 using the 
spectra with med(z) < 2.5; these are consistent to about 
la. Using the Q and QGS UVBs, the [O/C] values in- 
ferred from the redshift subsamples are marginally con- 
sistent with each other and with the full sample. 

4.2. tqvi versus thi for the full sample 

While the tovi/tciv ratios give the most direct con- 
straints on [O/C], it is also useful to examine tovi/thi, 
since comparing the simulated to the observed tovi/thi 
ratios gives an additional (but related) estimate of [O/C] 
(recall that our simulation reproduces the observed 

7"Civ( T Hl))" 

Figure 8 shows logTovi/ T Hi versus thi for our com- 
bined sample, in bins of redshift. Lines again connect the 
simulation points (with an overall scaling to best match 
the observations), which reproduce the observed trends 
in z and thi- The scalings correspond to best-fit [O/C] of 
0.69 ± 0.06, 0.19 ± 0.06, and 2.25 ± 0.06, for QG, Q, and 
QGS, respectively. For QG and QGS, the inferred [O/C] 
are consistent with the results found using tovi/ t civ; for 
our hardest UVB model, Q, the [O/C] values are more 
discrepant, but still within 1.5c. 

However, while the simulations reproduce the overall 
trends present in the data, there are some possible dis- 
crepancies. Although the logrni < 0.3 points are up- 
per limits, the simulations also appear to fall marginally 
above the data in this regime. At logrni ^ 1, on the 
other hand, the simulations slightly but significantly un- 
derpredict tovi- This can be seen more clearly in Fig- 
ure 9, where log tovi / Thi versus z in bins of thi is shown, 
exhibiting a clear discrepancy for points at high-z and 
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TABLE 2 
Best-fit Abundances 



UVB model 


[0/C] a 


x 2 /d.o.f. 


[0/C] b 


[Si/C] 


X 2 /d-o.f. 


[Si/O] 


[C/H] 

a 


=a+/3(z - 3)+7(log<5 - 

P 7 


0.5) 
X 2 /d-o.f. 


QG 
Q 

QGS 


0.66 ± 0.06 
0.06 ± 0.06 
2.23 ± 0.06 


42.7/82 
59.9/82 
26.7/82 


0.56 ±0.08 
0.06 ±0.08 
2.11 ±0.09 


0.77 ±0.05 

, 40+O.O5 
J--48-0.06 
26+ - 06 


65.7/115 
65.6/115 
73.8/115 


0.11 ±0.08 
1.42 ±0.08 
-1.97±0.11 


o 47 +0.07 

-2.91±0.07 
_4 14+0.06 
*' 14 -0.05 


nna+ 09 6 c i+ ' 10 
-0.06±0.09 0.17±0.08 
0-54±°;J° 1.31±0.07 


114.1/184 
113.8/184 
114.2/184 



Note. — The best-fit abundances and corresponding x 2 /dof from this work, and Papers II ([C/H]) and III ([Si/C]). Errors were computed by 
bootstrap-resampling the quasars used in the x 2 minimization. 



a Computed using tovi/tciv vs. Tcivl see § 4.1. 

^Computed using tovi/thi vs. thi for gas of density 6 < 10; see § 4.2. 




Fig. 8. — Rebinned median logTovi/Tn vs. log thi in bins of 
z for the combined QSO sample. The top two panels show bins 
centered at z = 2 and 3 with width Az = 1; the bottom panel shows 
combined data for all redshifts. Data points are plotted with la 
and 2<r error bars. The lines represent corresponding simulation 
points (with errors suppressed and with [O/C] chosen to minimize 
the x 2 ) using different UVB models, as indicated in the legend in 
the top panel. 



high thi- Indeed, if we consider subsamples above and 
below log thi = 1, we find that for logrui <1 we obtain 
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Fig. 9. — Rebinned median log(rovi/' r Hl) vs. z in cuts of thi 
for the combined QSO sample. Data points are plotted with 1 and 
2cr error bars, where green lines denote lower error of -oo. The top 
two panels show bins centered at log thi = —0.50 and 0.50 with 
width 1.0 dex; the bottom-left panel corresponds to a bin centered 
at logTciv = 1-75 with width 1.50 dex; the bottom-right panel 
shows data for all tciv combined. The lines are the same as in 
Fig. 8 



[O/C] values consistent those obtained from rovi A"civ 
(0.62 ± 0.07, 0.13 ± 0.07 and 2.18 ± 0.07 for QG, Q, 
and QGS respectively), but for log thi > 1 we obtain 
[O/C] of 0.92 ±0.10, 0.34±0.10, and 2.52±0.11 for QG, 
Q, and QGS respectively. This discrepancy in [O/C] for 
between low- and high-Tfji subsamples is significant at 
* 2.5 a for the QG and QGS models, and at 1.7 a for the 
Q model. 

Given the tight (but z-dependent) relation between thi 
and gas density (see Paper II, Fig. 2, and the upper 
axis in the top two panels of Fig. 8), it is useful also to 
divide our sample into high- and low-density subsamples. 
We have done this by recomputing [O/C] from rovi/ T Hi 
using only bins with thi corresponding to 6 < 10 or S > 
10. For S <10, we obtain [O/C] of 0.56 ± 0.08, 0.06 ± 
0.08, and 2.11 ± 0.09 (all consistent at la with the full 
sample), whereas for for S > 10 we obtain [O/C] of 1.02 
± 0.09, 0.51 ± 0.09, and 2.60 ± 0.09 for QG, Q, and 
QGS, respectively (all discrepant at « 4<r). 

This high-<5 and (less significantly) high-rui difference 
might have several causes. First, it might correspond to 
a genuine change in [O/C] with gas density. In this case, 
however, such an effect would also be expected in the 
trends of tovi/tciv vs. tqiv ~ aud no such effect is evi- 
dent. Thus, we consider it more likely that there exists a 
significant portion of hot gas - not present in the simu- 



Oxygen in the IGM 



9 




0.6 



Fig. 10. — Rebinned median log rovi/ r Silv vs. logrsnv in bins of 
z for the combined QSO sample. The two panels show bins centered 
at z = 2 and 3 with width Az = 1. Data points are plotted with 
1 and 2<r error bars. The lines represent corresponding simulation 
points (with errors suppressed and with [O/C] chosen to minimize 
the x 2 ) using different UVB models. The logrsuv = —1.35 Q data 
point is set to the minimum allowed tovIj logrovi=-6. 



lations - that contains O VI, but lacks H I and C IV. As 
suggested in §3.5 and Fig. 3, and further discussed below 
in §5.4, this might indicate the presence of a significant 
amount of collisionally ionized (T > 10 5 K) gas in the 
IGM. 

4.3. tovi versus Tsiiv for the full sample 

A third check on our measured [O/C] is provided by 
the rovi/ T Siiv ratio. Fig. 10 shows Tovi/ T Siiv versus 
Tsnv for our full sample, for two cuts in z. For each 
UVB the constant Si/C ratio derived in Paper III and 
the C distribution measured in Paper II are imposed on 
the simulations and the [O/C] ratio is varied to mini- 
mize the x 2 difference between the observed and simu- 
lated data points. Though the detection of rovi/T"Siiv 
is weak, the simulations appear to adequately represent 
the observations in the redshift range where rsnv is best 
detected, 2.5 < z < 3.5. From these fits we infer, for 
this redshift interval, [O/C] of 0.77 ± 0.19, 0.24 ± 0.20, 
and 2.17 ± 0.21 for UVB models QG, Q, and QGS, re- 
spectively; all are consistent with the results obtained 
from tovi/tciVj though all are somewhat higher. Be- 
cause Si IV probes higher density gas than C IV (sec 
Paper III), this again suggests that the simulations un- 
derpredict the amount of O VI in and near dense gas. 

4.4. [O/H] versus 8 from tovi/thi vs. thi 

As a final method, we can apply the "inversion 
method" developed in Paper II, to convert tovi/ t hi vs - 
thi into [O/H] vs. 6 by applying a density-dependent 
ionization correction. Then, using the measured distri- 
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Fig. 11. — Median oxygen abundance vs. overdensity S for 2=1.5- 
3.5. The inversion results for each QSO are binned in density bins 
of 0.125 dex. Data points are plotted with Icr error bars, where 
green lines denote lower error of -co. The solid line indicates the 
least-squares fit to the individual data points for S < 10. Dotted 
curves represent the lcr confidence limits, which were computed by 
bootstrap resampling the QSOs. The dashed line is the measured 
[O/H] utilizing the derived [O/C] results obtained from tovi/ t CIV 
and assuming the [C/H] distribution from Paper II. 



bution of carbon from Paper II, an independent measure- 
ment of [O/C] can be obtained. 

In Figure 11 we show the derived [O/H] versus 5 for 
our preferred UVB model QG, with data from all z com- 
bined. The data points from the individual quasar spec- 
tra (an example is shown in Fig. 4) have been binned in 
density bins of 0.25 dex. The solid line shows the least- 
squares fit to the data points with 8 < 10 and the dotted 
curves indicate the 1 a confidence limits, with the result- 
ing fit given in the upper left corner. The errors on the 
fits were determined by bootstrap resampling the QSOs. 
The dashed line is the value of [O/H] given by the de- 
rived [O/C] =0.64 result using tovi/ t civ, and assuming 
a [C/H] distribution from Paper II. 

For 5 < 10, the [O/H] derived from the ionization cor- 
rection agrees very well with that determined from both 
1 "ovi/' r civ, and tovi/ t hi- This strengthens the result 
from §4.1 and §4.2 that [O/C] is indeed constant for - 
0.5< log 6 < 1. 

However, for 5 > 10 the ionization correction results 
in substantially more [O/H] than that predicted using 
the previously-determined [C/H] distribution and a fixed 
[O/C] ratio. This is similar to the breakdown between 
the simulations and observed tovi/thi seen in Fig. 8, 
and the erroneously high [O /H] recovered for high-5 (see 
Fig. 5) when applying the ionization correction to the 
simulated spectra. Once again, this suggests the presence 
of collisionally ionized O VI. 

4.5. Systematic uncertainties 

The primary source of systematic uncertainty in, e.g. 
[O/C] or [O/H], is the complex modeling that must be 
performed to extract these values from the pixel correla- 
tions. The greatest combination of importance and un- 
certainty is clearly the uncertainty in the shape of the 
UVB; but this is discussed at length below in Sec. 5.1 so 
we here focus on other aspects of the modeling. 

The good agreement between the four methods we have 
employed indicates that the method is sound; but it is 
also clear that there are real differences between the uni- 
verse and our simulation. In particular, it is clear that 
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the real universe has an extra component of O VI at high 
density, almost certainly due to collisionally ionized gas 
that is not captured by those simulations. Nonetheless, 
if we exclude those high-Tni regions, the small discrepan- 
cies between [O/C] as measured using the different meth- 
ods indicate that such effects probably do not contribute 
more than ~ 0.1 dex uncertainty to our basic results. 

Another source of error that may be inaccurately as- 
sessed by our bootstrap resampling technique is that 
from continuum fitting. As shown in Table 1, our es- 
timated rms error in the O VI absorption region is 
~ 1 — 4%. To test the effect of this error on our re- 
sults, we have imposed an additional error on each ob- 
served spectrum on scales of 20, 80, and 320 A, for a total 
added rms error of 2%, then recomputed our results. We 
find that our best fit [O/C] from rovi versus tqw, and 
from rovi versus thi, are both within ~ 0.03 dex. Fur- 
ther, both the individual binned points and linear fits of 
the [O/H] values computed from the full sample (Fig. 11) 
are all affected by this continuum error to a lesser degree 
than the quoted random errors. Thus, we conclude that 
the continuum fitting error is not a significant systematic 
affecting our results. 

A final possible source of systematic error is that the 
O VI and C IV recombination rates used in the ver- 
sion of CLOUDY we have employed are too high by 
~ 50 — 75% compared with recent experimental values, 
for the temperatures relevant to low-density photoion- 
ized gas. (Savin, private communication). For the range 
of densities 1 < 5 < 40 we cover, this would imply a 
(density-dependent) correction of — 0.1 dex. This might 
change our overall fits by an amount comparable to the 
statistical errors, but is still much smaller than uncer- 
tainties stemming from the UVB shape, and cannot ac- 
count for the 0.5 dex of excess in tqvi at high thi > 10 
shown in Fig. 11. Combining these possible sources, we 
estimate probable systematic errors of ~ 0.2 dex in our 
basic [O/C] and [O/H] values; this uncertainty may be 
somewhat greater in subsets of the data, particularly at 
high-density. 

4.6. Densities and volume/mass fractions to which the 
results apply 

It is important to emphasize that each quoted result is 
sensitive to, and applies to, only a certain range of gas 
densities. At the upper end, our results nominally con- 
cern gas of up to 8 ~ 100, though (as noted at length 
above) the high-rni range of our data is likely to be af- 
fected by collisional ionization. That range is, however, 
not dominant: we have checked that if all pixels with 
thi > 30 are excluded from the analysis, the results 
given for [O/C] in columns 2 and 4 of Table 2 change 
only within the quoted errors. 

The lower end of the gas density range probed is most 
straightforward in results from tovi/thi (Figs. 8 and 11, 
which formally show OVI detections at la using all QSOs 
at z=2.5 — 3.5 for logTm > —0.1, or log 5 > 0.1 and 
confident detections at log 5 > 0.4. In Paper III, the Si 
abundance results were sensitive to logrni ^0.2 (log 5 > 
0.2), and in Paper II, C abundances were measured in 
much lower density gas. Thus our results provide indirect 
constraints on [O/C] and [O/Si] down to log 5 ~ 0.1 — 0.2. 

However, there are important caveats. First, the 
quoted results pertain to the the full density range 



probed and thus arc not necessarily very sensitive to 
the lowest densities. Second, because the pixel method 
only works if the element on the x-axis is more eas- 
ily detectable than the element on the y-axis, direct 
measurements of [O/C] and [O/Si] from tovi/tciv and 
T ovi/' r Siiv are dominated by much higher-density that 
gives Tciv^Siiv 3> —2.0; thus while these results are 
consistent with our indirect constraints, they do not ad- 
dress the (somewhat implausible) possibility that at low 
densities O, Si, and C come from completely different gas 
phases. On the other hand, at high densities it is quite 
possible that OVI and SilV emission are dominated by 
different phases, so the indirectly inferred [Si/O] values 
are probably both reliable and well-measured only in the 
moderate density range S ~ 5 — 10. 

Although mass or volume filling-factors corresponding 
to these results are not well constrained (see, e.g., Schaye 
& Aguirre 2005), the forewarned reader can convert the 
density range 5 > 2 correspond to into a volume using 
Paper II. 

5. ANALYSIS AND DISCUSSION OF RESULTS 

5.1. Relative abundances and the spectral shape of the 

UVB 

The best fitting metallicities and corresponding 
X 2 /d.o.f. from Papers II, III and this work are shown 
in Table 2 for each UVB model. Two interesting results 
stand out. First, for all UVB models carbon is under- 
abundant relative to both silicon and oxygen (being only 
marginally consistent with solar for UVB model Q). Sec- 
ond, all abundance ratios are sensitive to the UVB shape. 
A harder UVB results in a lower inferred [O/C] but a 
higher inferred [Si/C], making the [Si/O] ratio particu- 
larly sensitive to the spectral hardness of the UVB. 

The extreme sensitivity of the inferred [Si/O] ratio 
to the spectral shape of the UV background makes it 
possible to constrain feasible UVB models by making 
only weak assumptions about the [Si/O] ratio. Since Si 
and O are both a elements, they are expected to trace 
each other relatively well. For example, using the nucle- 
osynthetic yields of Portinari et al. (1998) and Marigo 
(2001) and a Chabrier (2003) initial mass function from 
0.1 — 100 M©, the [Si/O] ratio of the ejecta of a pop- 
ulation of age t > 10 s yr is predicted to be about 0.12 
and -0.03 for stars of solar and 1 percent solar metal- 
licity, respectively. This agrees well with observations 
of metal-poor stars, which find [Si/0]« (Cayrel et al. 
2004). 

Tallying the results of this work with those of Papers 
II and III yields [Si/0]= 0.11 ± 0.08, 1.42 ± 0.08, and 
-1.97 ±0.1 for UVB models QG, Q, and QGS, respec- 
tively. Thus, our preferred model, QG, is nicely consis- 
tent with the expectations, but models Q and QG lead 
to inferred [Si/O] ratios that are highly inconsistent with 
both nucleosynthetic yields and observations of metal- 
poor stars. (Assuming [Si/C] ~ 0.5 and [Si/O] ~ in 
the Q background, for example, raises the \ 2 of the fits 
in Figures 6 and 8 by 65 and 92, respectively; requiring 
this for QGS likewise raises the x 2 by 138 and 81.) We 
conclude that the UVB has a spectral shape similar to 
that of model QG. 

While our result using the QG UVB are broadly con- 
sistent with the abundance ratios in metal-poor stars and 
in yield calculations, the [O/C] and [Si/C] may be some- 
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what high, by ~ 0.1 — 0.3 dcx. This is comparable to our 
systematic errors, but nevertheless interesting if taken 
seriously. 

For example, the models of Nomoto et al. (2006) that 
include the contributions of hypernovae (defined as su- 
pernovae with kinetic energy > 10 x that of normal core- 
collapse SNe) produce [O/C] w 0.6, and [Si/C] « 0.65, 
in agreement with our results. 

5.2. Implications for cosmic abundances 

In Paper II we combined the median [C/H](<5, z) with 
the width cr([C/H])(<5, z) of the lognormal probability dis- 
tribution of [C/H] for —0.5 < log 5 < 2.0 to determine 
the mean C abundance versus 5. This was then inte- 
grated over the mass- weighted probability distribution S 
(obtained from our hydrodynamical simulation) to com- 
pute the contribution by gas in this density range to the 
overall mean cosmic [C/H]. Assuming that [O/C] is con- 
stant over this density range 9 we obtain, for our fiducial 
UVB model QG, [0/H]= -2.14 ±0.14, corresponding to 

n OJGM ~ 3.3 x io-6 10 [o/h ]+2 .i (Jh_"j (5) 

Extrapolating our [C/H] and [O/C] results even further 
to the full density range of the simulation would yield 
values w 0.2 dex higher but with more uncertainty, as we 
have argued that our results are unreliable at the highest 
densities. 

Note that these results are relatively insensitive to the 
UVB (unlike those for fisuGM m Paper III) because for 
a harder UVB, the inferred [C/H] increases, while [O/C] 
decreases. For our quasar-only model Q, these effects 
almost entirely cancel, yielding [0/H]w —2.3, and an 
^O.igm value 30% lower than for model QG. Note also 
that these estimates include the oxygen that resides in 
gas that is observable in C IV and O VI, but they do 
not include oxygen in intergalactic gas that is very hot 
(T > 10 5 K) or very cold (T < 10 4 K) and shielded 
from ionizing radiation. 

If, following Bouche et al. (2007), we take Oz.igm = 
^O,igm/0.6, we then infer an intergalactic metal reser- 
voir of Oz,igm ~ (4.3 — 6.8) x 10~ 6 . This can be compared 
to their estimate of the total z « 2 — 3 "metal budget" of 
flz ~ 2 - 3 x 10~ 5 , indicating that ~ 15 - 35 % of metals 
produced prior to z — 3 reside in the component of the 
IGM that is studied here. 

5.3. Previous Searches for Oxygen in the IGM 

Previous studies have explored oxygen abundances in 
the IGM using both line-fitting (e.g., Rauch et al. 1997; 
Dave et al. 1998; Carswell et al. 2002; Bergeron et al. 
2002; Simcoe et al. 2004; Bergeron & Herbert-Fort 2005) 
and pixel optical depths (Dave et al. 1998; Schaye et 
al. 2000a; Telfer et al. 2002; Aracil et al. 2004; Pieri & 
Haehnelt 2004). 

Previous pixel studies did not attempt to convert their 
O VI detections into oxygen abundances; but we can 
compare to their recovered optical depths. Both Schaye 
et al. (2000a) and Aracil et al. (2004) claim detection 
of O VI down to thi ~ 0.2. Using our combined data 

9 Note that this is an extrapolation beyond the range 2 < 8 < 10 
over which we have reliably measured Oxygen abundances. 



set binned in density (see Fig. 11), we obtain l-er de- 
tections down to about the mean density (thi ~ 0.5 at 
z = 2.5). While O VI is in principle an excellent tracer of 
metal in very low-density gas, in practice we find that for 
the higher redshifts where low densities are more easily 
probed, H I contamination is severe. Thus we are not 
in practice able to constrain metals in underdense gas 
as claimed in previous studies, in spite of a large sample 
and improved techniques of removing contaminants. On 
the high- thi side, both Aracil et al. (2004) and Pieri & 
Haehnelt (2004) exclude pixels saturated in H I, so can- 
not probe log thi ^ 0.5. (This accounts, for example, 
for our detection of O VI in Q1422+230, while Pieri & 
Haehnelt (2004) had no detection in the same QSO). The 
study of Schaye et al. (2000a) did probe high thi, where 
their results are broadly consistent with ours. 

The studies of Carswell et al. (2002), Simcoe et al. 
(2004), Telfer et al. (2002), and Bergeron & Herbert- 
Fort (2005) did perform ionization corrections and we 
can compare our abundance determinations (relatively 
directly) to theirs. Carswell et al. (2002) assumed to 
abundance of oxygen relative to carbon to be solar and 
inferred metallicities for various UVB models and a num- 
ber of O VI absorbers in TVhi ~ 10 15 cm~ 2 systems at 
z w 2. They found that for relatively hard UV back- 
grounds, comparable to our model Q, the ionization mod- 
els yielded densities in agreement with theoretical pre- 
dictions for self-gravitating clouds with the observed H I 
column densities (tih ~ 10~ 4 cm -3 ; Schaye (2001)) and 
metallicities of 10~ 3 — 10~ 2 solar, in excellent agreement 
with our measurement of [O/C] « [C/H] « —2.6 for model 
Q at z = 2 and log 5 = 1.5. We note that if Carswell et 
al. (2002) would have allowed oxygen to be overabundant 
relative to carbon, they would have found that softer 
UVB models are required to obtain density estimates 
that agree with theoretical expectations for gravitation- 
ally confined clouds. 

Telfer et al. (2002) employed the Faint Object Spec- 
trograph on HST to search for O V in QSO spectra from 
redshift 1.6< z < 2.9. The O VI/O V ratio found by 
the survey favored a UVB background similar to our Q, 
which they use to derive a metallicity of —2.2 <[0/H] 
< —1.3. While somewhat higher than our value, much 
of the difference may be attributable to their use of the 
78th percentile in the O V/H I ratio. Assuming that 
[O/H] has a scatter at fixed density similar to that in 
[C/H] (see Paper (II)), this would correspond to a me- 
dian of - 0.5 dex less, or -2.7 <[0/H] < -1.8, in fairly 
good agreement with our numbers. 

Bergeron & Herbert-Fort (2005) divide their sample 
into "metal-poor absorbers" with N(0 VI) /N (H I) < 
0.25, which they take to be predominantly photoionized, 
and "metal-rich absorbers" with N(0 VI)/7V(H I) > 
0.25, for which they assume a hotter phase. For the 
metal-poor systems they use a "hard" UVB and assume 
[0/C]= to derive a range of -3.0 < [O/H] < -1.0, 
and for the metal-rich phase they infer a median [O/H] 
« —0.80 to —0.33, depending upon the assumptions re- 
garding the ionization balance. Combining their samples 
they estimate cosmic density fio.iGM ~ (1.6—4.4) x 10~ 6 , 
corresponding to —2.4 < [O/H] < —2.0 if divided by the 
cosmic gas density = 0.045. While precise compar- 
ison is difficult, these numbers are consistent with our 
corresponding estimates of [0/H]rs —2.3, or fiojGM ~ 
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3.3 x 1CT 6 using the Q UVB. 

Simcoe et al. (2004) assume [O/C]=0.5 and 0.0 for 
UVB backgrounds comparable to our QG and Q, for 
2-2 < z < 2.8, so (as for the above studies) compar- 
ing derived [O / C] values is less useful than several other 
points of comparison. First, in both backgrounds, their 
dependence of [O/H] upon S is similar to that found for 
[C/H] vs. S in Paper II, consistent with a constant [O/C] 
value. Second, Simcoe et al. (2004) find that there is 
a clear jump in in the median [O/H] at S ~10, while 
a corresponding jump is not seen in [C/H], similar to 
our results in §4.4. They also interpret this as possibly 
indicating that stronger absorbers are physically more 
complex or multiphased. Third, Simcoe et al. (2004) 
compute an overall contribution fio.iGM ~ 1-0 x 10~ 6 , 
using their "hard" background; this would correspond 
to a cosmic average contribution of [O/H] =—2.6; these 
numbers are w 0.3 dex lower than our values, but this 
should be regarded as good agreement given the number 
of assumptions made in each computation. 

Finally, the studies of Rauch et al. (1997)Dave et al. 
(1998) are similar to ours in employing simulated spectra 
to attempt to match observed O VI absorption and thus 
constrain [O/C] and [O/H]. Rauch et al. (1997) gener- 
ated simulated spectra from a constant-metallicity sim- 
ulation and compared ionic ratios to extant data, infer- 
ring [C/H] « —2.5 and evidence for overabundance of 
Si and O relative to carbon. Dave et al. (1998) used 
Q 1422+230 and a quasar-only Haardt & Madau UVB 
model like 'Q'. The found that their data is consistent 
with [O/C] w and [C/H]« —2.5; these are quite con- 
sistent with our results using the Q UVB. As an alterna- 
tive interpretation, they note that a softer UVB would 
give high [O/C] (more characteristic of Type II super- 
nova yields), but also lower [C/H], again consistent with 
our findings; however they interpret this softness as due 
to patchy reionization, whereas we favor its explanation 
by the contribution of galaxies to the UVB. 

5.4. Collisionally Ionized Gas 

As discussed above, the difference between the inferred 
[O/H] for 5 > 10 and 5 < 10 is probably due to collision- 
ally ionized gas, for reasons sketched in § 3.5. Fig. 3 
shows that for T 3> 10 5 K, collisional ionization domi- 
nates and the optical depth ratios become independent 
of the density. For T <C 10 5 K, on the other hand, the 
O VI/C IV and the O VI/H I ratios both drop rapidly 
with increasing density. Consequently, these ratios can 
be many orders of magnitude higher in hot, dense gas 
than in warm, dense gas. The lower the density, the 
smaller the differences become. In fact, at high densities 
the O VI fraction in warm gas is too small for O VI to be 
observable. Therefore, any O VI at the redshift of very 
strong H I absorption is likely to arise in a different phase 
than the associated H I and possibly even the associated 
C IV. The O VI phase must either have a much lower 
density or a much higher temperature. 

Because the fraction of hot gas in our simulation is 
small at all densities, we effectively assume the gas to be 
photo-ionized when we compare with synthetic spectra 
and when we correct for ionization (as in Fig. 11). In the 
latter case, we also implicitly assume that O VI and H I 
absorption arise in the same gas phase. 

Hence, our results from § 4.4 suggest the existence of 



a detectable amount of O VI enriched hot (T > 10 5 K) 
gas associated with strong H I absorption. A possible 
physical explanation is that such O VI systems coincide 
with outer regions of high z galactic halos, where the 
effects of galactic winds may dominate the heating pro- 
cess. If the temperature exceeds 10 5 K in these regions, 
then this gas would contain significant O VI, while lack- 
ing C IV and H I. The latter two would then arise in a 
cooler gas phase, which has to be fairly dense in order 
to account for the strong H I absorption. The high den- 
sity of the cooler phase implies that it would not produce 
significant photo-ionized O VI. 

These results and inferences are consistent with other 
observational studies of O VI. Carswell et al. (2002) and 
Bergeron et al. (2002) both find that the majority of 
the detected absorption systems had temperatures de- 
termined from the line widths too low for collisional ion- 
ization, but cannot rule out higher temperatures for some 
absorbers. Indeed, the study of Simcoe et al. (2002) find 
their detected high-column density O VI lines associated 
with strong H I absorbers are broad enough to be consis- 
tent with collisional ionization. As noted above, Simcoe 
et al. (2004) find a jump in [O/H] at 6 ~10 interpretable 
as a transition to a regime in which collisionally ionized 
gas affects the abundance inferences. 

While collisionally ionized gas complicates oxygen 
abundance inferences, the flip-side is that O VI then pro- 
vides an important probe of hot, enriched IG gas that is 
difficult to detect in C IV (e.g., Paper II). Indeed, hy- 
drodynamical simulations by Theuns et al. (2002) and 
Oppenheimer & Dave (2006) predict that a significant 
portion of C IV is collisionally ionized. Hence, it would 
thus be very interesting to compare such simulations, 
employing O VI as a metallicity tracer, with the obser- 
vations analyzed here. 

6. CONCLUSIONS 

We have studied the relative abundance of oxygen in 
the IGM by analyzing O VI, C IV, Si IV, and H I pixel 
optical depths derived from a set of high-quality VLT 
and Keck spectra of 17 QSOs at 2.1 < z < 3.6, and 
we have compared them to realistic, synthetic spectra 
drawn from a hydrodynamical simulation to which met- 
als have been added. Our fiducial model employs the 
ionizing background model ("QG") taken from Haardt 
& Madau (2001) for quasars and galaxies (rescaled to re- 
produce the observed mean Lya absorption). The simu- 
lation assumes a silicon abundance as calculated in Paper 
III, [Si/C]= 0.77 ± 0.05, and a carbon abundance as de- 
rived in Paper II: at a given overdensity 5 and redshift z, 
[C/H] has a lognormal probability distribution centered 
on -3.47 + 0.65(log(5 - 0.5) and of width 0.70 dex. The 
main conclusions from this analysis are as follows: 

• For 1.9< z <3.6, r H i < 10, and 5 <10 (when 
smoothed on the scale of the H I absorption, 10 — 
10 2 kpc), the fiducial simulation utilized in Papers 
I, II, and III consistently agrees with the observed 
tovi/tciv(tciv), tovi/thi(t H i), and (to a lesser 
degree) rovi/Tsav^siiv)- Fitting t vi/Vciv(tciv) 
yields a constant [O/C] = 0.66 ± 0.06, with es- 
timated systematic errors within ±0.2 dex. Con- 
verting the observed tovi/thi(thi) into [0/C](<5) 
using the ionization correction method of Paper II 
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further supports these results. 

• The relative abundances [O/C] and (especially) 
[O/Si] are sensitive to the UVB shape. We find 
that our fiducial (Haardt & Madau 2001, quasars 
and galaxies) spectrum gives reasonable results for 
both, but that significantly softer or harder UVBs, 
such as the Haardt & Madau 2001 quasar-only 
UVB, give results that arc highly inconsistent with 
both theoretical yields and observed abundance 
ratios in other low-metallicity environments, and 
should not be considered tenable. 

• Our results, both from applying the ionization cor- 
rection and from comparing the simulations to the 
observations, suggest no evolution in [O/H] over 
the rcdshift range 1.9< z <3.6, but a strong de- 
pendence on 5. Both results are consistent with 
those found in Paper II for [C /H] . 

• For thi > 10 and 8 > 10 the value of [O/C] (derived 
by comparison to the simulations) is inconsistent 
with that found at lower densities, and ~ 0.5 dex 
higher than that predicted using the carbon distri- 
bution of Paper II and a density-independent [O/C] 
value. This might in principle suggest a density- 
dependent [O / C] ratio, but we favor the interpreta- 



tion that a fraction of the high-J O VI absorbing gas 
is collisionally ionized, and that this leads to an er- 
roneously large ionization correction in this regime. 
This interpretation is supported by our simulated 
spectra as well as by the observation that O VI 
lines associated with strong H I absorbers tend to 
be broader than those associated with weak H I 
systems (Carswell et al. 2002; Simcoe et al. 2002; 
Bergeron et al. 2002). 
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